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Abstract: 

This communication is an enquiry into the circumstances under which concurrence 
and phase entropy methods can give an answer to the question of quantum entanglement 
in the composite state when the photonic band gap is exhibited by the presence of pho- 
tonic crystals in a three-level system. An analytic approach is proposed for any three-level 
system in the presence of photonic band gap. Using this analytic solution, we conclusively 
calculate the concurrence and phase entropy, focusing particularly on the entanglement 
phenomena. Specifically, we use concurrence as a measure of entanglement for dipole emit- 
ters situated in the thin slab region between two semi-infinite one-dimensionally periodic 
photonic crystals, a situation reminiscent of planar cavity laser structures. One feature 
of the regime considered here is that closed-form evaluation of the time evolution may be 
carried out in the presence of the detuning and the photonic band gap, which provides 
insight into the difference in the nature of the concurrence function for atom-field cou- 
pling, mode frequency and different cavity parameters. We demonstrate how fluctuations 
in the phase and number entropies effected by the presence of the photonic-band-gap. 
The outcomes are illustrated with numerical simulations applied to GaAs. Finally, we 
relate the obtained results to instances of any three-level system for which the entangle- 
ment cost can be calculated. Potential experimental observations in solid-state systems 
are discussed and found to be promising. 

PACS numbers: 42.50.Dv, 03.65.Ud, 03.67.Mn 



1 Introduction 

A structure in which the dielectric constant varies periodically is called a photonic crystal. 
One of the most interesting properties of a photonic crystal is the existence of a photonic 
band gap PQI2] • Radiation with a frequency that lies within the band gap cannot propagate 
in the photonic crystal structure. Photonic crystals are usually viewed as an optical analog 
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of semiconductors that modify the properties of light similar to a microscopic atomic 
lattice that creates a semiconductor band-gap for electrons [HEIII]. Photonic band gap 
crystals offer unique ways to tailor light and the propagation of electromagnetic waves and 
have caused growing interest in recent years because it offers the possibility of controlling 
and manipulating light within a given frequency range through photonic band gap pH 0] . 
Photonic band-gap materials have attracted much attention in recent years for theoretical 
and practical importance in fundamental science and application H] . The atom-photon 
interaction in photonic band gap materials [4_ has been found to exhibit many interesting 
new phenomena such as photon-atom bound states [Sj, spectral splitting [Hj, quantum 
interference dark line effect [7j, phase control of spontaneous emission [S], transparency 
near band edge P, and single-atom switching [TUj . 

In a parallel development, considerable work was done recently on entanglement prop- 
erties jllj . The detection of entanglement is one of the fundamental problems in quantum 
information theory. From a theoretical point of view one can try to answer the question 
whether a given entirely known state is entangled or not, but despite a lot of progress in 
the last years [11,12], no general solution of this problem is known. In experiments, one 
aims at detecting entanglement without knowing the state completely. Bell inequalities 
and entanglement witnesses ^3] are the main tools to tackle this task. Interestingly, 
the concurrence of the ground state which is related to the entanglement of formation, 
has been shown to be strongly affected at the critical point [16J. More precisely, in the 
one-dimension, it has been shown that the derivative of the concurrence with respect to 
the coupling constant diverges at the transition point, although the concurrence itself is 
not maximum. These pioneering results raise the question of the universality of these 
behaviors. Actually, the lack of exact solutions especially in higher dimensions implies 
a numerical treatment which often restrict the study to a small number of degrees of 
freedom. 

Heisenberg's uncertainty relations had tremendous impact in the field of quantum 
optics particularly in the context of the construction of coherent states and also for differ- 
ent physical systems as well as the reconstruction of quantum states. The minimization 
problem of finding the number-phase uncertainty state has been considered and mini- 
mum uncertainty state relations between number and phase uncertainty are presented 
[T7] . Many authors argued that ^2]; the Heisenberg inequality is too weak for practi- 
cal purposes, which led them to the establishment of information theoretic uncertainty 
relations. 

Our aim of the present paper is to consider the dynamics of a system of three-level 
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atoms with dipole interaction in presence of the photonic band gap and study the concur- 
rence and the entropic uncertainty relation for number and phase. With applying some 
approximations, one can deal with the quantization of the electromagnetic field modes 
of a homogeneous, but anisotropic medium which can then be made to form one of the 
sandwich layers in the slab structure under consideration involving two semi-infinite peri- 
odic photonic crystals. With the electromagnetic modes quantized, one can evaluate the 
entanglement degree, and explore its variations with the controllable parameters of the 
system. To reach our goal we have to find an exact analytic solution of the time dependent 
Schrodinger equation of the system. We show that a reasonable amount of entanglement 
can be achieved in a system of three-level atoms with dipole interaction in presence of 
the photonic band gap and essentially we establish deeper connections between entropic 
uncertainty relations and entanglement. 

The organization of this paper is as follows: in section 2, we give an overview of 
effective medium approach and dispersions, followed by subsection 2.1 where we introduce 
our Hamiltonian model and give exact analytic solution for the Schrodinger equation in 
the frame of the dressed state formalism. In section 3, we employ the analytical results 
obtained in section 2 to investigate the properties of the entanglement degree due to the 
concurrence, and classify the behavior in several parameter regimes assuming that the 
electromagnetic field is in a coherent state in subsection 3.1. In section 4, we essentially 
establish deeper connections between entropic uncertainty relations and entanglement. 
Numerical results for the phase entropy are discussed in the subsection 4.1 for two different 
cases; one is the resonant and the other is the off-resonant case. The prospects for 
experimental observation of our predictions are analyzed in section 5. Finally, a summary 
of the main points of this work ends the paper and a few avenues for further investigations 
are indicated in section 6. 

2 Effective medium approach 

The effective-medium approach can be applied to situations in which all three regions of 
the structure possess frequency-dependent dielectric functions. In fact, the rapid pace of 
the technological progress in solid-state quantum computing gives one a hope that the 
specific prescriptions towards building robust qubits and their assemblies discussed in this 
work can be implemented in future devices. In this regard, very promising fields where 
the concept of nonlinear localized modes may find practical applications is the quantum 
computation of photonic band gap materials, periodic dielectric structures that produce 
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many of the same phenomena for photons as the crystalline atomic potential does for elec- 
trons [4,18]. Nonlinear photonic crystals (or photonic crystals with embedded nonlinear 
impurities) create an ideal environment for the generation and observation of nonlinear 
localized photonic modes. Much theoretical work has been done on the properties of fi- 
nite one dimensional photonic band gap (PBG) crystals including recent calculations 
of the thermal emissivity of such one- dimensional structures The strong angular 

dependence of the gap effect with a one-dimensional structure has motivated successful 
experimental work with three-dimensional structures [2U|. in particular, the existence 
of such modes for the frequencies in the photonic band gaps has been predicted j2l] for 
2D and 3D photonic crystals with Kerr nonlinearity. Nonlinear localized modes can also 
be excited at nonlinear interfaces with quadratic nonlinearity [22J, or along dielectric 
waveguide structures possessing a nonlinear Kerr-type response [23] • 

Z + 




Figure 1: A schematic representation of the dielectric slab between two photonic crystals 
occupying the regions z < and z > r. 

The system that we consider here consists of a dielectric cavity occupying the region 
< z < r and the photonic crystals occupy the regions z > r and z < 0. For long 
wavelength fields and in the effective medium approach, the photonic crystal has the 
optical characteristic of a uniaxial medium [3J. Moreover we shall specialize to uniaxial 
media, so that our system has only two principal axes with the z— axis as the optical axis. 
In this case, the components of the dielectric tensor appropriate to the photonic crystals 
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where, e = e e", e z = e e z . The dielectric tensor components for the two semi- infinite crys- 
tals can be written in the following forms = (r] 1 di+ri 2 d2)/di2, € z i = 7 7i J 72 G ^i2/ (^7i<^2 + 
r) 2 di), and ej| = (773^3 + 774^4) / 'd 3Al e z2 = V3V4 d 3A/(r] 3 d 4 + r] 4 d 3 ), where the d {j = di + dj, 
the subscripts 1 and 2 on ej and e Z i refer to the first and second photonic crystal. The 
dielectric functions r] 1 and r] 2 , one or both of which may be frequency dependent. The 
photonic crystals are treated using the effective medium approach, which pertains to any 
layer structure formed by alternate periodic stacking of two types of layers of locally 
isotropic materials of thicknesses di and d 2 (see figure 1). 

In this paper we are concerned with the interface polaritons which are characterized 
by imaginary wave vectors normal to the interfaces such that the waves are decaying with 
distance from the interfaces at z = and z = r into the outer regions and are hyperbolic 
in the slab [3J 1221 1231 ■ To see the salient features of the effective medium description we 
shall ignore retardation effects, which amounts to ignoring throughout (u)/c) terms. In 
this case dispersion relation for the surface polaritons takes the form 

u ♦ u ( k ° ( fc i/ e i) + ( fc 2/e 2 ') \ 

k s r = arctanft x — - — — — f . (2) 

The dispersion relations obtained from the Maxwell wave equation of this system lead to 
two distinct equations k 2 s = — u> 2 e s /c 2 , kf = e\ k?Je z i — uo 2 ^} /c 2 , where s refers to 
the slab cavity. 

It is important to note that infinite and semi-infinite photonic crystals have the same 
band structure [26 . The only difference is the existence of surface modes in the case 
of semi-infinite structure. The main feature of all ID photonic crystals is that although 
forbidden gaps exist for most given values of the tangential component of the wave vector 
(k), there is not an absolute nor complete photonic band gap if all possible values of the 
tangential component of the wave vector are considered 4J. Having determined the modes 
we can now quantize the fields associated with these modes using the usual quantization 
procedure |3] the single-mode quantized field takes the form 

E(x,t) = Eoa(k\\) exp[i(k\\.x - urt)] + H.C., (3) 



5 



where E is the strength of the electric field, k\\ is the wave vector, x is the position 
operator and a the annihilation operator. 

2.1 The model and methods of solution 

Accurate potentials are of course required for a quantitatively correct prediction of the 
behavior and properties of real quantum systems. However, even qualitative conclusions 
drawn from simulations employing inaccurate or invalidated potentials can be problematic. 
The most appropriate form of the potential depends largely upon the properties of interest 
to the simulators. Now we consider the interaction of the abovementioned modes with a 
three-level atom in three different configurations, namely, V— , Lambda- and cascade-type. 
The transition in the 3-level atom is characterized by the dipole matrix element A^. The 
operator Su describes the atomic population of level \i)a with energy ujj, (j = a, b, c) and 
the operator Sij, (i ^ j) describes the transition from level \i)a to level \])a- The total 




Figure 2: The V-type, A-type and S-type three-level atom interacting with a bimodal field. 
The levels \1)a, |2)a, and |3)a, have the energy values huii, hix>2 and hu3, respectively. The 
transitions \1)a — > \2)a , and \2)a — ► \3)a, are coupled to two intra-cavity different 
modes Si and &2 with eigenfrequencies Qi and 1^2- The detunings of the levels \1)a, |2)a, 
and (2)^4, |3)a, are Ai = U\ — uj 2 — ^i and A 2 = u-s — uj 2 — ^2, for V-type, and S-type, 
while Ai = oj 2 — u)\ — fii and A 2 = 0J2 — ^3 — ^2 for A-type. 

Hamiltonian of this system is H = H Q + Hint- The 3 eigenstates, |£J and corresponding 
eigenenergies, ati are assumed to be known. The total wave-function may be expanded in 
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terms of the known eigenstates, namely 



\m) = Mt)\ti) + Mt)\t*) + Mt)\t*)- (4) 



With atomic units, using Schrodinger equation, we obtain the coupled equations for our 
three-level system, namely 

dt 



i ?Ml =rjAj (t)+ y £ i H jh Ak(t), (5) 



fe=i 



where H = r i an d Hjk = Hint |Cfc) • These equations are exact for any three- 
level atom. In the interaction picture, let us consider a three-level system described, in 
an appropriate rotating frame, by the Hamiltonian 

Hint = Ai5ll + A 2 S , 3 3 + A2I-R1S21 + A32-R2532 + X21R1S12 + Xl 2 RlS 23 . (6) 

The atom-field couplings A y - are given by = Y/jl^.E, where E is the quantized electric 
field given by equation (3) and ■ is the matrix dipole moment coupling between the state 
i and j. The Y factor accounts for local field effects and is given by Y = 3e s (u;) / (2e a {u) + 
1), where e s (u) is given in equation (1). It is easy to write Ay in the following form 

= 3e,(u;) {uj/ujt) 2 - (u L /u T ) 2 
13 2e.{w) + l' (lu/cut) 2 -r/ 2 ' 1 ' 

where rj 2 = [ 2e s (uj) (ul/ut) 2 + l]/[2e 8 (u) + 1]. The transitions between the three lev- 
els may occur in three different configurations depending upon the relationship between 
the energies Ei,E 2 and E 3 of levels 1,2 and 3. The possible configurations are [SJJ 
(i) the V-type corresponding to E 2 < E\ < E 3 , (ii) the A— type or Raman configura- 
tion corresponding to E\ < E 3 < E 2 and (iii) the S— type or ladder-type corresponding 
E\ < E 2 < E 3 . Each of the two pairs of levels can be coupled by only one-mode or 
two-mode. The field operators in the abovementioned three types are (i) F\ — a\ F 2 — b 
for V-type, (ii) F\ = a, F 2 = & for A -type and (iii) F% = a, F 2 = b for S-type with a = b 
if both pairs of levels are coupled by the same mode. 
In order to solve equations (5), we assume that [2*7] 



G(t) =A(t)+xB(t)+yC{t), (8) 



which means that 



% — 7T- = r i + ViV) < A(t) + ■ — —B[t) + ■ — —C(t) } , 9 

dt I n + v*y r x + v\y 1 
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where V\ and v 2 are given using equations (5) and (6). We seek G{t) such that iG(t) = 
zG(t). This hold if 

v 2 x + r 3 y r 2 x + v 2 y t 

y = — ; — — x = — ; — — > z = n+v x y. 

ri + v{y r\ + v{y 

After some algebra this leads to a cubic equation which has three eigenvalues Xi{y,i) 
which determine the Z{. There are also three corresponding eigenfunctions Gj(t) = 
Gj(0) exp(— izjt), where 

Gj(t) = MjxA(t) + M j2 B(t) + M j3 C(t), (10) 

where 

/ 1 xi yi^ 

1 x 2 y 2 . (11) 

\ 1 ^3 V3 J 

Now, we express the unperturbed state amplitude A(t),B(t) and C(t) in terms of the 
dressed state amplitude Rj 

3 3 

Fi(t) = M ^Gj(t) = V G i(°) exp(-^t), (12) 
j=i j=i 

-^1,2,3 {t) = A, B, C. Using the above equations, we can write 

A(t) = [(x 2 y 3 - V2X 3 )e- iZlt + {x 3 y t - y 3 xi)e- iZ2t + {x x y 2 - yix 2 )e~ iZ3t ] , 
B{t) = i [(y 2 - y 3 )e~^ + (y 3 - y^e^ + ( Vl - y^e'^] , (13) 

C(t) = -^[(x 2 -X3)e- l2lt -(x3-x 1 )e- t22 '-(x 1 -x 2 )e-^ t ], 

where D = det(M) = x\y 2 + x 2 y 3 + x 3 y\ — x\y 3 — x 2 y\ — x 3 y 2 . We have thus completely 
determined the dynamics of a three-level system in the presence of photonic crystal. 

The picture in this case is of the three-level system in the presence of photonic band 
gap and the detuning, rather than the usual picture of the three-level Jaynes-Cummings 
model (JCM) system. The important point to note here is that, using the above an- 
alytic approach, any three-level Hamiltonian is likewise exactly solvable, with precisely 
similar eigenvectors and eigenvalues that are obtained directly using equations (4) and 
(6). In Ref. [2H1 a n analytic approach is proposed for three-level systems, based on the 
Riccati nonlinear differential equation. However, the solution obtained is valid only in 
certain situations. On the other hand, our analytic approach removed the restriction that 
considered in the previous work and this solution is valid for any three-level system. 

Next, we discuss a frequently encountered phenomena of particular interest in which 
we define the entanglement measure of the present system. 
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3 Concurrence 



Quantum entanglement has recently been attracted much attention as a potential resource 
for communication and information processing [2H] . Entanglement is usually arise from 
quantum correlations between separated subsystems which can not be created by local 
actions on each subsystem. The concept of concurrence originates from the seminal work 
of Hill and Wootters JH] where the exact expression of the entanglement of formation of 
a system of two qubits was derived. They showed that the entanglement of formation, 
an entropic entanglement monotone, is a convex monotonic increasing function of the 
concurrence. 

It has been shown that the concurrence of a mixed two-qubit state, C(p AB ), can be 
expressed in terms of the minimum average pure-state concurrence, C {\iPab))-> where 
the minimum is taken over all possible ensemble decompositions of Pab- So that, the 
concurrence is defined of a mixed state p for 2x2 quantum systems, in the following form 

m 

C(p) = max (a i - cr 2 - er 3 - er 4 ) , (14) 

where the <7; are the square roots of the eigenvalues of the product matrix Q, the singular 
values (by convention sorted in descending fashion), all of which are non-negative real 
quantities 

Q = ^(Jy ^ Oy^fp, (15) 

Oy is the well-known Pauli matrix, and yfp is any matrix satisfying ^fp = yfp* . The 
importance of this measure follows from the direct connection between concurrence and 
entanglement of formation Ef 

E f (p) = -p + \np + -p_\np_, (16) 

where 

^ = \ (l ± yfi ~ C{ P f) . (17) 

One can prove that p is separable if and only if the concurrence is zero. 

Let us now turn our attention to the definition of the concurrence of a pure state [SHj 
on a (N x K) —dimensional Hilbert space 3ft = ® 31k- The flip operator F acting on 
an arbitrary Hermitian operator A on 9ft can be written as 

F(A) :=A + (trA)U - (tr N A) ®U K -U N (tr K A), (18) 

where tr N and tr k the partail traces over $l N and 3?^-, respectively. We denote by IIw 
and Uk the identity on and 3?^-, respectively. The expectation value F(p^) , 
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where = is non-negative for all pure states and equals zero if and only if 

is a product state. This allows to define the concurrence of any arbitrary bipartite pure 
state as fUJl 

cm) = y/m nv> m 

= ^2{(m 2 -tr(p%)), (19) 

where p N = trx (pw,) is the reduced density operator of dimension N. For a normal- 
ized state, (if)\if)) = 1, it interpolates monotonously between zero for product states and 
^ 2(N-i) max j ma ]]y entangled states. 

To investigate the concurrence for the system under consideration, we have to evaluate 
the reduce atomic density matrix p A = tr F p(t), which can be written as 

p A = E pJ0<»1+ £ Pij \i)(j\, (20) 

i=l,2,3 i,j=l,2,3,v£j 

where p^t) = (i\p A (t)\j), i,j — 1,2 and 3. Using equations (19) and (20), we can write 
the concurrence in the following form 



C(m = j2 £ (puPh-P^P*)- ( 21 ) 

V i,j=l,2,3,i^j 

Although the concurrence and therefore the results we obtain are not restricted to the 
standard one-mode three-level system, we will use that language throughout most of the 
paper. 

Having specified the various photonic crystal and field amplitude parameters, we will 
present in the following subsection the results of our numerical analysis of the concurrence. 



3.1 Numerical results 

For applications in real systems, we consider the dipole emitters with frequencies in the 
reststrahl band of GaAs. In this subsection we will discuss the time dependence of 
the concurrence, which considered as an entanglement measure. We will consider the 
commonly used state as initial condition for the cavity field: the coherent state, which 
may be applicable in different situations. As might be expected, the behavior of the 
three-level system changes dramatically depending on the initial field state. Throughout 
this subsection the quantity to be examined is the concurrence C . 

In figure 3, we present the oscillatory behavior of the concurrence C against 
the scaled time Ait and the mean photon number n for different values of the detuning 
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Figure 3: The evolution of the concurrence C (\ip)) as a function of the scaled time \±t 
and the mean photon number n. The parameters are eo = 10.89, rj = 1.085, uj/ujt = 2, 
uj /uj t = 1, Hu L = 36.29meV, Hu T = 33.25meV, d x = 500 A, d 2 = 300A, e x = 9, e 2 = 1.3, 
c?3 = 500A, <i 4 = 400A, €3 = 10, 64 = 1.5 and L = 1.5c?, and different values of the 
detuning parameter, where A = for Fig. 3a and A = 5Ai for Fig 3b. 



11 



parameter, where A = for Fig. 3a and A = 5Ai for Fig 3b. We consider a specific system 
in which the cavity is taken as GaAs with e = 10.89, rj = 1.085, uo/uot = 2, Uq/ujt = 1, 
Hlu l = 36.29meV, ftw T = 33.25meV. The photonic crystals parameters are given by the 
arbitrary set, di = 500A, d 2 = 300A, e x = 9, e 2 = 1.3, d 3 = 500A, d A = 400A, e 3 = 10, 
64 = 1.5 and L = 1.5c?. The general behavior due to the coherent state of the field does not 
contain any surprises it is quite broad, corresponding to the standard quantum limit. The 
value of concurrence at the first maximum is 1, which is quit remarkable, see figure 3a. 
After the time goes on, we see that the maximum value of the concurrence decreases with 
small amplitude of the oscillations. As the mean photon number increased, the number 
of oscillations decreased. 

The effect of the parameter which describes the mismatch between the atomic fre- 
quency and the mean frequency of the cavity mode has been considered in figure 3b. We 
set the other parameters as the same as in figure 3a, and A = 5Ai. As A is increased 
the behavior of the three-level system becomes increasingly erratic. Shorter revival times 
cause successive revivals to overlap and interfere so that the time evolution appears ir- 
regular. The detuning parameter at which irregularity emerges is closely tied to the 
mean-photon number: the higher the mean-photon number, the smaller the detuning 
needed to produce irregular behavior. Larger detuning also results in decreased revival 
amplitude due to the larger number of frequencies in the sum, which causes the rephasing 
to be less complete. However, a signature of the revivals persists as a return to the bare 
Rabi frequency even at mean-photon number high enough that the behavior looks ran- 
dom and the revival amplitude is essentially washed out. From our further calculations 
(which are not presented here), we point out that as we increase the value of the detuning 
one can see that the revival time is also prolonged, however the period of fluctuations is 
decreasing. Detuning affects the revival time by elongating it and the maximum value 
of the entanglement degree becomes smaller and smaller. Similar to the case of a two- 
level atom, detuning shifted the atomic occupation probability around which it oscillates 
upward meaning that the energy is stored in the atomic system. 

Now we will turn our attention to the effect on the concurrence of the mode frequency 
uj/ujt- In particular, we consider eo = 10.89, i] = 1.085, uoq/lot = 1 and for different 
values of the scaled time, where Ait = tt/2 for Fig. 4a and Ait = 37r/2 for Fig 4b. Our 
particular observation is the maximum entanglement occurs near the band edges, which 
corresponds to u — 1.085ujt- Near the band edges the wave vector parallel to the interface 
reaches its maximum value, and this corresponds to the first two relatively small peaks 
around the point 1.085. In the gap region or the reststrahl region of the GaAs system 
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Figure 4: The evolution of the concurrence C as a function of the mode-frequency 
uj/ujt and the mean photon number n for different values of the scaled time, where 
\ x t = tt/2 for Fig. 4a and A a t = 3vr/2 for Fig 4b. 
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no electromagnetic fields can propagate and coupling is therefore suppressed. The extra 
peaks around the point 1.085 are attributed to local field effects and can be understood 
from looking at equation (7) where \j has a pole at 77 = uj/uot- One has to bear in mind 
that the above calculation did not take into explicit account the spatial dependence of 
the coupling parameters. Therefore, a more careful calculation would have to take into 
account the nonstationary property of the present system. The above model calculations 
suggest that physical parameters such as mode frequency, mode-atom coupling and cavity 
dielectric have important effects on the entanglement. One can see that the oscillations 
collapse after few Rabi periods and after an interval of time in which the concurrence is 
constant, the oscillations reappear again. This revival then collapses and a new revival 
begins. 

This behavior highlights once again the role of the functional form of the modified 
Rabi frequencies in controlling the time evolution of the concurrence. Rabi frequencies 
which obtained in the present model are similar to that obtained from the standard three- 
level model but involving a frequency-dependent dielectric function. An important point 
to keep in mind when comparing the results presented here with results from the usual 
three- level system in the absence of the photonic band gap is that: they are give a different 
feature relative to the entanglement. This raises an interesting question: can one use the 
present system in building quantum logic gates? Calculations and detailed discussion of 
this issue will be presented in a forthcoming paper. 



4 Phase entropy 

One of the most striking features of quantum mechanics is the property that certain 
observable cannot simultaneously be assigned arbitrarily precise values. This property 
does not compromise claims of completeness for the theory, since it may consistently be 
asserted that such observable cannot simultaneously be measured to an arbitrary accuracy 
[TT] . The Shannon entropies associated with the photon number distribution P m and phase 
probability distribution P(9,t), 

P m = (m\p(t)\m), 
P(e,t) = (B\p(t)\B), (22) 
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where \m) is the Fock state and \6) is the phase state, are given respectively by [T7j 

oo 

Rn = — Pm hi Pm, 

= - J (P(6, t) In P(9,t))d6. (23) 

The entropic uncertainty relations for the number and phase distribution determine the 
lower bound on the sum of the Shannon entropies Rn and R^ : 

Rn + R^p > ln(27r). (24) 

This equality is satisfied by a Fock state for which R^ = and R^ = ln(27r). Other 
physical states give an entropic sum greater than ln(27r). Specifically, for a coherent state 
we find that the sum is ln(e7r) for the mean photon number greater than one, i.e. 

Rn + R^ > ln(e7r). (25) 

The lower bound for the position-momentum entropic uncertainty relation is also given 
by right-hand side of this equation i.e ln(e7r). 

The single-mode of the Pegg-Barnett phase formalism which of interest in the field of 
quantum optics can be constructed from the single-mode phases [2] to take the form 

P(6,t) = lim ( S -±l)(9 m \ p (t)\9 rn ), (26) 



\8 m ) is a phase state of the mode, 

\°rn) = ^=Y j e me -\n) 1 (27) 

where m = 9 + ^y, and m = 0, l,...s, and 9 arbitrary. Equation (26) defines a 
particular basis set of (s + 1) mutually orthogonal phase states. 

Using the standard procedure [HI] , the phase probability distribution, the expectation 
value and the variance of the Hermitian phase operator may be obtained for the field. 
Since the coherent field at t = belongs to a class of partial phase states, we have chosen 
the reference phase 6*o as #o = (3 — jv^, and introduced the new phase labels ( = m — |s 
where m = 0, 1, 2, s. Then as s tends to infinity the summation may be transformed 
into an integral after replacing by 6, and by dO. This leads to continuous phase 
probability distribution, where 

P(0,t) = -L ( 1 + 2 jr {A n , m (t) cos[9(n - m)\ + B n>m (t) sm[9(n - m)}} J , (28) 

\ n>m / 
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Figure 5: P(9,t) against 9 and the scaled time Ait. The parameters are eo = 10.89, 
rj = 1.085, u/ut = 2, ujo/ujt = 1, hjJL = 36.29meV, Hut = 33.25meV, d± = 500A, 
d 2 = 300A, e x = 9, e 2 = 1.3, d 3 = 500A, d A = 400A, e 3 = 10, e 4 = 1.5 and L = 1.5c/, where 
(a) A = and (b) A = 5A X . 
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where A n>m (t) and B n>m (t) are given by 




Re {A n (t)A* m (t) + B n (t)B* m (t) + C n (t)C* m (t)} , 
Im {A n (t)A* m (t) + B n (t)B* m (t) + C n (t)C* m (t)} . 



(29) 



The phase probability distribution is normalized according to j^^P(9 ) t)d9 = 1. 



4.1 Numerical results 

In what follows we shall display some general arguments based on the equality sign in 
the Heisenberg uncertainty relations that to demonstrate the phase entropy of a general 
three-level system in the presence of photonic band gab when the initial state of the field 
is assumed to be in a coherent state. 

In figure 5a, we have plotted the phase probability distribution P{9,i) as a function 
of the scaled time Ait and 9 taking into consideration the presence of the photonic band 
gap. For example at time Ait = we realize that the phase distribution P(9 ) t) starts 
with a single-peaked structure at 9 = corresponding to the initial coherent state. Then 
as the time develops the peak splits into two peaks moving into two opposite directions. 
However the amplitudes of the split peaks fluctuate in time giving a top like shape until 
the two peaks reach the values 9 = ±n at middle of the revival time but in this range 
the amplitudes of the peaks do not show any fluctuations. The picture changes greatly 
as time develops further (say A x t > 40) where we find that the two-peak profile breaks 
up into multi peak with reduction of the amplitudes of these peaks. Thus the phase 
distribution shows diffusion as well as bifurcation. Different features are visible when we 
consider the off-resonant case and the behavior of the phase probability distribution is 
changed dramatically (see figure 5b). In this case we observe that there is a diffusion of 
the peaks at earlier time. 

In figure 6, we consider the behavior of the phase probability distribution against the 
mode frequency uo/oot and 9 for the same parameters as in figure 5, while in this figure, we 
keep the the scaled time Ait fixed, where, Ait = tt/2 for figure 6a and Ait = 3ir/2 for figure 
6b. One may clearly see that the phase probability distribution is discontinuous near the 
band edges. This corresponds to the zero value at the point 1.085. We can prove, in an 
analogous manner to the equation (7), that the uj/ujt = 1.085 is a pole of the atom- field 
coupling can not avoided. It is interesting to see that, the phase probability distribution 
does not depend on the mode frequency for a fixed value of 9, except at the point 1.085. 
As the time increased, the only difference is that, the phase probability distribution peak 
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Figure 6: P(9,t) against the mode-frequency uj/uot and 9. The parameters are the same 
as in figure 5, but (a) Ait = n/2 and (b) Ait = 3n/2. 
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splits into two peaks moving into two opposite directions, keeping the symmetry around 
the point 9 = 0. 
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Figure 7: The phase entropy (a) and the number entropy Rn (b) as functions of the 
scaled time Ait with an initial coherent state of the radiation field with n = 20 based on 
the exact numerical results due to equation (23). We consider the same parameters as in 
figure 1. 

In figure 7, we plot the number entropy Rn and the phase entropy as functions of 
the scaled time Ait. The initial state of the field is considered as a coherent state. We 
specifically present the results for the same values of figure 5. It should be noted that at a 
special choice of the mean-photon number parameter, the situation becomes interesting, 
where the Rabi frequency has a minimum value at n. In this case we find that the general 
behavior of the entropies Rjsr and R^ and with an initially coherent field exhibit irregular 
structures instead of the regular structure resembling those manifested by the number 
or vacuum states cases. Here it is interesting to note that the periodic oscillations are 
observed for a short period of the interaction time only. When we consider smaller mean- 
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Figure 8: The phase entropy (a) and the number entropy Rn (b) as functions of the 
scaled time X\t where, n = 10 and the same parameters as in figure 1. 
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photon number, the regularity behavior of the oscillations in the entropies Rn and 
are still obvious (see figure 8) where we have considered the initial mean photon number 
n = 10. However, the number of oscillations is increased. Also it is interesting to point 
out that at the revival time optimal phase entropy is attained in all the cases which means 
that the atom has achieved an almost pure state, this has been observed all through our 
figures. The number and phase entropic uncertainties for a weak coherent state follow 
those of underlying states superposition. 
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Figure 9: The phase entropy (a) and the number entropy Rn (b) as functions of the 
mode-frequency u/ujt, where, ft = 20, Ait = n/2. 

In figure 9 we plot the entropies R^ and against the mode-frequency u in units 
of ujt for different values of the scaled time. Now, where the atom-field coupling is 
proportional to \j this explains the origin of the second peak in this figure. It is interesting 
to note the dependence of these entropies on the mode-frequency, with different values of 
the scaled time. We wonder, as a possible generalization of this concept, whether there 
exist another family of similar oscillations if we consider two-qubit system. In such a 
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case, the properties of these systems would probably be of interest, in order to bring 
further insight and knowledge about entanglement and quantum logic gates for multi- 
partite systems. As the scaled time increased, a characteristic feature of the entropies 
and is quit interesting, where more oscillations exist, also, only around the resistable 
region for the number entropy but less number of oscillations exist for the phase entropy 
(see figure 10). Far from the resistable region the entropies behavior observed here does 
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Figure 10: The phase entropy R^ (a) and the number entropy Rn (b) as functions of the 
mode-frequency uj/ujt, where, n = 20, = 37r/2. 

not depend on the mode frequency and the intensity of the initial field mode. 

Given a 3-D photonic crystal with a complete gap, one has the possibility of introducing 
a defect in the structure which will create a localized state in the gap. If this is a point-like 
defect then the photon mode will be completely localized about a point. In figure 10, we 
show the zero point associated with the defect created by removing a small amount of 
dielectric from one of the vertical dielectric columns of the crystal structure. The resulting 
defect mode has a state near mid gap. One feature that should be highlighted in this 
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context is the appearance of a frequency gap between the pair of interface dispersion. This 
gap is present only when the two photonic crystal regions are different, and disappears 
when they are identical. 

5 Experimental prospects 

The perfect semiconductor crystal is quite elegant and beautiful, but it becomes ever 
more useful when it is doped. Likewise, the perfect photonic crystal can become of even 
greater value when a defect is introduced [32J. The point to make about photonic crys- 
tals is that they are very empty structures, consisting of about 78 empty space. But 
in a sense they are much emptier than that. They are emptier and quieter than even 
the vacuum, since they contain not even zero-point fluctuations within the forbidden fre- 
quency band. Our model system consist of a three-level atom located inside a photonic 
band gap material. There are several ways of placing such an atom inside a photonic 
crystal. From a material standpoint, it is possible to dope an existing photonic band gap 
material using ion beam implantation methods. For instance, it has recently been shown 
that Er 3+ ions implanted into bulk silicon exhibit sharp free-atom- like spectra |3*2^ IHH] . 
Intense temperature-dependent photoluminescent (PL) at 1.54 fim is observed in the sys- 
tem at low temperatures (when the host material is crystalline, Er-related PL is quenched 
at temperatures above 80K so that it cannot be detected at room temperatures). This 
wavelength is particularly significant because it corresponds to the minimum absorption 
of silica fibre-based optical communication system. Because the PL at 1.54/xm is due to 
the spin-orbit split 4 ii3/2 — > 4 I15/2 of 4/ electrons in the Er 3+ ions which are shielded 
by outer 5s 2 5p 6 shells, the influence of the host lattice on the luminescence wavelength 
is weak. (The key to the success of erbium is that the upper level of the amplifying 
transition 4 /i3/2 is separated by a large energy gap from the next-lowest level 4 /i5/2 so 
that its lifetime is very long and mostly radiative. In spite of the screening of the atomic 
transition by the outer shells, it is likely that thermal phonons in the silicon host would 
cause significant dephasing of the quantum degrees of freedom within the erbium 4f shell. 
Consequently, such a system must be cooled to liquid helium temperatures. Such experi- 
ments appear to be nearly within the reach of current technology. Although it has not yet 
been demonstrated, the system consisting of a multi-level system coupled to a multi-mode 
appears to be another potential candidate for achieving new features. Such systems are 
potentially interesting for their ability to process information in a novel way and might 
find application in models of quantum logic gates. Therefore, atoms or trapped ions + 
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cavities in a presence of photonic band gap represent, in our opinion, a very promising 
system for quantum information processing. 

6 Conclusion 

In this communication the quantum electrodynamic properties of a three-level atom em- 
bedded in a photonic band gap material were investigated. We have focused on the ap- 
plication of the effective-medium theory to the present problem in a nanoscale dielectric 
cavity QED situation. The effective-medium approach can in fact be applied to situa- 
tions in which all three regions of the structure possess frequency-dependent dielectric 
functions. Specifically, the combined effects of coherent control by an external driving 
field and photon localization facilitated by a photonic band gap on entanglement from a 
three-level atom embedded in a photonic band gap material were examined. Exact solu- 
tions of the wave function in the Schrodinger picture have been obtained within rotating 
wave approximation. In particular, we have chosen to focus on three-level system coupled 
to a single mode. Observation of the three-level system may offer some insight into the 
quantum nature of the resonator, just as atoms provide a sensitive probe for the nonclas- 
sical nature of electromagnetic fields. The observation of revivals, which are a strictly 
nonclassical phenomenon, would give evidence for the quantum nature of the quantum 
system. 

The results point to a number of interesting features, which arise from the variation 
of the adjustable parameters of the system, namely, the mode-frequency, dipole vector 
orientation, dipole position within the slab, the slab width, and the photonic crystal pa- 
rameters: layer widths and dielectric functions. Our investigations for the entanglement, 
collapse-revival phenomena, and phase and number entropic uncertainty relations in the 
presence of the photonic band gap as compared with the usual three-level model are 
summarized as follows:- 

i) The concurrence behavior is reflect the pattern of collapse and revival which is 
qualitatively similar to that of the usual three-level model but with reduced amplitude. 
In case of a smaller mean photon number and for initially excited atom the usual pattern 
in the three-level model of collapse and revival changes to rapid fluctuations of interference 
patterns for all time considered. In this way, our concurrence function contains all the 
information necessary to identify the entanglement of a given state. Nevertheless, it 
depends on the particular choice of the mode-frequency. 

ii) The phase entropy can be used to measure entanglement of the system presented 
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here with explicitly atom-field coupling in the presence of photonic band gap. We would 
like to point out that the phase Shannon entropic considered for the presented model has 
not been treated in this manner before. 

iii) The photonic band gap introduces sudden changes in the concurrence and phase 
entropy due to the variation of these quantities with mode frequency. This feature at- 
tributed to the fact that in the photonic band gap region electromagnetic modes are not 
allowed to propagate into the dielectric slab and hence no interaction can take place in 
this region. Theory predicts analytically this behavior for a GaAs system at u — tjut- 

Finally, we emphasize the fact that without any conditions it was possible to obtain 
exact analytic solution which reproduce the most important features of the three-level 
atom interacting with a cavity one- or two-mode in the presence of photonic band gap. 
A similar set of equations have been derived in [2H] for a three-level system using some 
approximations, based on the Riccati nonlinear differential equation. In contrast, the 
method used here gives exact analytic solutions without any conditions. 
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